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ABSTRACT 

Several  algorithms  and  techniques  have  been 
developed  to  estimate  the  horizontal  directionality 
of  ambient  noise  from  horizontal  line  arrays.  In 
most  cases,  the  two-dimensional  (2-D) 
representation  of  the  array  beam  pattern  is 
adequate  to  estimate  the  ambient  noise  horizontal 
directionality.  However,  when  significant 
amounts  ofthe  noise  are  distributed  out  of  the 
horizontal  plane,  a  2-D  representation  of  the  noise 
field  can  be  misleading.  An  algorithm  is 
described  herein  that  utilizes  the  3-D  conical 
nature  of  the  beam  response  patterns  of  the  towed 
array  and  an  iterative  technique  to  estimate  the  3- 
D  ambient  noise  directionality.  The  algorithm  is 
described  and  results  from  actual  at-sea 
measurements  are  presented.  This  work  was 
supported  by  funding  from  the  ASW 
Environmental  Acoustic  Support  (AEAS) 
Program  of  the  Office  of  Naval  Research. 

INTRODUCTION 

An  algorithm  to  estimate  the  three  dimensional 
(3-D)  ambient  noise  field  has  been  developed  that 
utilizes  horizontal  line  array  beam  noise  data  on 
multiple  array  headings  as  the  basic  data  input. 
Several  measurement/processing  techniques  have 
been  developed  to  estimate  the  horizontal 
^directionality  of  ambient  noise  frorr,  horizontal 
tine  arrays1-2.  In  most  cases,  the  2-D 
representation  of  the  array  beam  pattern  is 
adequate  to  accurately  estimate  the  horizontal 
directionality1.  When  significant  amounts  of  the 
noise  are  distributed  at  elevation  angles  well  off 
the  zero  elevation  angle,  use  of  the  3-D  structure 
of  the  beam  patterns  and  a  priori  information 
about  the  noise  vertical  arrival  structure  greatly 
improve  the  estimate  of  the  horizontal 
directionality  (2-D  noise  field)2.  Unfortunately, 
the  3-D  noise  field  is  not  estimated.  When  the 
vertical  component  of  the  noise  changes  as  a 
function  of  azimuth,  as  is  generally  the  case  when 
shipping  noise  is  an  important  contribution,  the 
true  3-D  conical  beam  response  pattern  needs  to 


be  implemented  so  that  the  noise  energy  will  be 
accurately  distributed  in  elevation  as  well  as  in 
azimuth.  The  algorithm  described  herein  is  an 
iterative  process  that  uses  the  true  3-D  conical 
nature  of  the  horizontal  line  array  beam  response 
pattern  during  the  3-D  ambient  noise  directionality 
estimation  process. 

In  order  to  formulate  the  ambient  noise  3-D 
estimation  algorithm,  the  primary  equation  that 
transforms  the  spherical  coordinates  of  the  noise 
field  into  the  conical  coordinates  of  the  array  will 
be  presented.  The  formulation  of  the  technique  is 
provided  in  a  sufficiently  general  nature  to  handle 
deviations  from  the  array  axis  including  an 
arbitrary  array  tilt.  Finally,  actual  results  from  a 
recent  noise  field  measurement  exercise  will  be 
used  to  demonstrate  the  capability  of  the  3-D 
ambient  noise  field  estimation  algorithm. 

BACKGROUND 

The  line  array,  because  of  its  relatively  low 
cost  and  ease  of  deployment,  is  often  used  to 
measure  the  undersea  ambient  noise. 
Furthermore,  the  towed  horizontal  line  array, 
with  its  added  advantage  of  mobility,  has 
achieved  a  high  degree  of  utilization  by 
researchers.  The  major  disadvantage  of  the  line 
array  is  that  it  can  only  distinguish  the  acoustic 
arrival  angle  relative  to  the  axis  of  the  array.  It 
cannot  uniquely  determine  the  spherical 
coordinates  of  the  arrival  angle  of  the  sound.  The 
resulting  ambiguity  in  angle  of  arrival  describes  a 
cone  with  apex  at  the  center  of  the  array  and  its 
center  axis  oriented  along  the  axis  of  the  array. 
The  arrival  angle  of  the  sound  relative  to  the  array 
axis  is  half  of  the  cone  apex  angle.  This  is 
illustrated  by  the  representation  of  the  towed  array 
in  Fig.  1  and  the  superimposed  3-D  beam 
response  patterns.  Because  of  the  vertical 
stratification  of  the  acoustic  propagation  field  in 
the  ocean,  the  various  beams  of  the  line  array 
sample  different  noise  field  and  acoustic 
propagation  regimes,  i.e.  SOFAR  channel 
(RRR)  and  boundary  interacting  (RSR). 
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Th<?  complex  nature  of  the  coupling  of  the 
reception  characteristics  of  the  line  array  (conical 
beams  with  the  spherical  structure  of  the  noise 
field  is  partially  illustrated  by  Fig.  1.  However, 
it  is  illustrated  better  in  Fig.  2  in  which. the  noise 
sphere  is  represented  as  a  flat  surface  and  the 
beam  response  "footprints"  Of  both  a  horizontal 
line  aiYay  on  a  heading  of  270°  with  a  few  degrees 
'  of  array  tilt  (solid  curves)  and  a  vertical  line  array 
(dashed  horizontal  lines)  have  been 
superimposed.  Only  the  footprints  of  the 
horizontal  line  array  have  both  elevation  and 
azimuttir  angle  dependencies.  Therefore,  it  should 
be  possible" to  estimate  both  the  azimuth  and  the 
elevatfon  angle  dependence  of  the  noise  from 
horizontal  line  array  data  when  measurements  are 
made  on  several  different  headings.  A  slight  tilt 
of  the  array  as  illustrated  in  Fig.  2  is  desirable  to 
eliminate  upward-downward  ambiguities.  It  is 
very  difficult  to  balance  the  array  well  enough  that 
it  always  tows  perfectly  horizontal.  Hence,  some 
degree  of  tilt  is  generally  a  natural  consequence  of 
the  measurement  process. 

The  undersea  ambient  noise  Field  at  a  given 
location  can  be  expressed  as  a  function  of  time 
and  spherical  angles  in  the  following  form: 

K(9.Q,t)  =  N(0,<t>)+^(©,<M)+e(0><M)  (1) 

where 

t  =  time 

0  =  azimuth  angle  in  spherical  coordinates 

0  =  elevation  angle  in  spherical  coordinates 

N(9,<?,t)  =  measured  noise  power  at  a  given 
frequency 

N(9,0)  =  pseudostationary  background 
ambient  noise 

p(0,o,t)  =  time  dependent  component  of  the 
ambient  noise  due  to  fluctuation  in  propagation, 
noise  source  movement,  changes  in  noise  source 
levels. 

E(0.o,t)  =  errors  in  measurements  due  to 
flow  noise,  non-  straight  array,  calibration  errors, 
etc. 

In  an  ideal  measurement  e(0,<J>,t)  is 
eliminated.  Various  measurement  and  data 
processing  techniques  have  been  developed  to 


minimize  ,s  effects  in  the  results.  p(0,<{>,t)  can 
be  ignored  if  only  the  pseudostationary 
component  N(©,<{>)  is  desired.  It  is  this 
pseudostationary  component  that  can  be  expected 
to  repeat  on  a  seasonal  basis  year  after  year,  when 
the  noise  source  distributions  and  acoustic 
propagation  conditions  are  again  the  same.  It  is 
this  term  that  can  be  used  in  the  predictions  of  the 
noise  for  future  times,  and  it  is  also  this  term 
which  the  technique  presented  herein  attempts  to 
estimate. 

The  technique  for  estimating  N(0,<())  is  based 
on  beam  noise  data  measured  by  a  horizontal  line 
array.  This  can  be  expressed  mathematically  as: 

T  2n  n/2 

Ri,j=4  Md0  J  N(0,d>,t) 

0  0  -rc/2 

bi(0-rj,<D)  cosO  d<D  (2) 

where 

Rj  j  =  beam  output  noise  powers  averaged 
over  a  time  T, 

i  =  beam  number, 
j  =  an:ay  heading  number, 

Tj  =  the  array  heading  angle,  and 
bj  (0,  <|>)  =  power  response  of  ith  beam. 

APPROACH 

The  horizontal  line  array  can  determine  the 
arrival  angle  X  of  a  signal  relative  to  the  axis  (or 
line)  of  the  array,  but  it  cannot  determine  the 
absolute  orientation  in  space.  Hence,  all 
orientations  give  the  same  response  as  long  as  the 
arrival  angle  relative  to  the  array  axis  is  the  same. 
The  resulting  beam  pattern  of  such  an  array 
response  is  conical,  which  intersects  a  concentric 
sphere  in  a  circle  as  illustrated  in  Fig.  3.  The 
noise  Field  on  the  other  hand,  is  more  naturally 
described  in  terms  of  spherical  coordinates.  The 
.'mathematical  representations  of  the  conical  beams 
of  the  line  array  and  the  spherical  nature  of  the 
noise  Field  are  basically  incompatible.  This 
incompatibility  must  First  be  overcome  before  an 
efficient  algorithm  can  be  developed  to  estimate 
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the  spherical  noise  N(0,<J))  from  conical  data  Rjj. 
'litis  incompatibility  is  overcome  by  Equation  (3) 
which  transforms  the  noise  power  at  spherical 

angles  0  and  <j>  into  noise  at  the  conical  angle  p: 

P  =  cos-1  (cos©  cosd?  cosX  +  sin  [cos'l(cos0 
cos<I>)]  sinX  sin  [tan"l(sin<b/cosd>  sin©)] }  (3) 

The  estimation  of  N(0,<(>)  is  an  iterative  one 
that  begins  with  assuming  an  estimated  noise  field 

N0(©,4>).  'Fife  simplest  possible  form  is  a 
constant,  whicjt  corresponds  to  an  isotropic  noise 
field  (e.g.  No(0,({))  =  C).  Next  N.-  (0,<p)  is 
convolved  with  the  beam  patterns  of  the  array, 
using  equation  (3),  to  get  estimates  R’j  j,  for  Rj  j. 
The  differences  Dy  =  Rj  ;  -  R’y  are  used  to 
either  add  noise  to  or  subtract  noise  from 
No(0,<t>)  at  the  appropriate  locations  0  and  $. 
When  this  has  been  done  for  all  i  beams  and  j 
array  headings,  a  new  estimate  Ni(0,(j>)  has  been 
created.  This  process  of  convolution  and  noise 
field  modification  is  continued  until  the 
differences  fall  below  a  predetermined  threshold, 

and  then  the  final  noise  field  estimate  Nf(0,4>)  is 
accepted  as  the  estimate  for  N(0,<t>). 

RESULTS 

Beam  noise  data  from  a  line  array  on  different 
headings  in  the  undersea  ambient  noise  field  were 
prpeessed  by  the  3-D  noise  field  estimation 
algorithm.  The  results  are  presented  in  Fig.  4. 

The  noise  levels  as  a  function  of  elevation  angle  <J> 

and  azimuth  angle  0  are  plotted  in  gray  scale 
contours  with  the  level  increasing  with  darkness. 

The  horizontal  directionality  of  the  noise  N(0)  is 
one  statistical  characterization  of  the  noise  field 
that  a  2-D  estimation  algorithm  would  attempt  to 
generate.  In  the  case  of  the  3-D  noise  field,  it  can 
be  obtained  by  merely  summing  over  all  vertical 

angles  <}>  at  each  azimuth  angle  0  (bottom  plot  in 
Fig.  4).  The  curves  in  the  rectangular  plot  to  the 
right  in  Fig.  4  are  vertical  profiles  of  the  total 
noise  and  of  two  small  azimuthal  sectors,  one  in  a 

direction  of  high  level  noise  (0  =  290  deg)  and 
one  in  a  direction  of  low  level  noise  (0  =  150 


deg).  The  vertical  profile  of  the  total  noise  is 
generally  obtained  from  measurements  by  a 
vertical  line  array.  In  this  case,  it  was  obtained 

by  summing  over  all  azimuth  angles  0  at  each 

elevation  angle  <{>.  The  other  two  vertical  profiles 

are  for  small  sectors  in  0  that  could  not  be 
measured  by  either  a  vertical  line  array  or  a 
horizontal  line  array.  In  fact,  the  measurement  of 
those  profiles  would  normally  require  a  high 
resolution  volumetric  array  with  an  extremely 
marrow  beam  pattern,  something  that  presently 
“does  not  exist. 

SUMMARY/CONCLUSIQNS 

A  technique  for  using  towed  horizontal  line 
array  beam  noise  data  measured  on  multiple 
headings  of  the  array  has  been  developed.  It 
utilizes  an  iterative  technique  that  compares 
measured  beam  noise  data  with  corresponding 
results  obtained  from  convolving  the  beam  noise 
patterns  with  an  estimate  of  the  3-D  noise  field  to 
guide  modifications  of  the  noise  field  estimate. 
Once  the  final  estimate  is  obtained,  it  is  then 
possible  to  investigate  characteristics  of  the  noise 
field  that  are  essentially  invisible  to  the  array  that 
acquired  the  original  beam  noise  data,  such  as  the 
vertical  directionality  of  the  noise  in  azimuthal 
sectors  that  otherwise  would  require  a  high 
resolution  volumetric  array  to  measure. 
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Figure  1.  Representation  of  a  towed  horizontal  line  array  measuring  undersea  ambient  noise. 
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Figure  2.  Representation  of  a  spherical  noise  field  N(0,  <|>)  and  the  conical  "footprints"  of  the 
vertical  line  array  (VLA)  and  the  tilled  horizontal  line  array  (HLA)  on  a  plane  surface. 
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Figure  3.  Coordinate  system  for  the  noise  sphere  and  the  conical  beams  of  the  line  array. 
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